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ABSTRACT 


A computer  program  is  described  for  calculating 
time-dependent^  laminar  natural  convection  flows  in 
a room.  The  complete  Fortran  language  program  as 
used  on  the  NBS  Univac  1108  computer  is  given.  A 
discussion  of  the  program  and  instructions  for  its 
use  are  facilitated  by  the  aid  of  a flowchart  and  an 
example. 
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COMPUTER  PROGRAM  TO  CALCULATE 
NATURAL  CONVECTION  FLOWS 
IN  A ROOM 


1.  INTRODUCTION 


This  report  provides  a description  of  the  computer  program  EXCONl 
and  instructions  for  its  use.  The  program  was  written  in  ASA  Fortran 
language  (reference  1)  and  coded  for  use  on  the  NBS  Univac  1108  computer. 
The  program  was  developed  to  calculate  time-dependent^  laminar  natural 
convection  flows  in  a room.  The  room  is  in  the  form  of  a vertical 
circular  cylinder^  with  the  heat  source  which  induces  the  flow  located 
in  the  center  of  the  floor  (see  figure  1) . The  computer  program  is 
based  upon  mathematical  principles  which  are  discussed  in  references  2 
and  3.  A comparison  of  several  programs  for  calculating  flows  is  presented 
in  reference  2,  which  highlights  the  particular  advantages,  as  well  as  the 
limitations,  of  EXCONl  (known  therein  as  numerical  method  V) . An  appli- 
cation of  EXCONl  is  presented  in  reference  3,  in  which  calculated  results 
are  compared  with  the  experimental  flows  of  reference  4. 

Figure  1 illustrates  a cylindrical  room  of  height  a and  radius  b. 

In  the  center  of  the  floor,  a disk  of  radius  c is  held  at  a temperature 
At  above  the  remaining  walls  of  the  enclosure.  The  flow  is  assumed  to  be 
axisymmetric  two-dimensional;  i.e.,  without  variations  in  the  azimuthal 
direction.  A grid  is  constructed  as  shown  of  M+-1  horizontal  lines  and 
N+1  vertical  lines  spaced  at  intervals  of  AX  and  AR.  The  program  EXCONl 
is  employed  to  calculate  the  temperature,  vorticity,  stream  function  and 
vertical  and  radial  velocities  at  the  intersections  of  these  grid  lines 
(called  grid  or  mesh  points) . By  rotation  about  the  centerline,  the 
resulting  fields  of  temperature,  velocity,  etc.,  are  translated  to  any 
azimuthal  plane. 


2.  DISCUSSION  OF  EXCON  1 


An  annotated  listing  of  the  complete  program  EXCONl  is  presented 
in  figure  3.  A flowchart  of  the  main  program  is  presented  in  figure  2. 
The  next  two  sections  describe  the  main  program  and  the  subroutines, 
respectively. 

2.1  Main  Program 

The  logical  sequence  of  the  main  program  (pp  1-7  of  figure  3)  can  be 
grasped  most  readily  by  first  considering  the  flowchart  in  figure  2.  A 
list  of  input  parameters  is  read  in  on  data  cards  (an  example  of  which  is 
shown  on  p 16  of  figure  3)  as  follows; 
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aspect  ratio  of  enclosure  = b/a. 
relative  size  of  heat  source  = c/a. 

Grashof  number  = gpATa^/v^^  where  g = gravity^  P = volume 
expansion  coefficient  of  fluid^  and  V = kinematic  viscosity 
of  fluid. 

Prandtl  number  of  fluid  in  enclosure  = v/K^  where  K.  = 
thermal  diffusivity. 

number  of  vertical  mesh  points  minus  one. 
number  of  radial  mesh  points  minus  one. 

Time  level  at  start  of  calculation. 

Time  step  for  first  time  advancement.  If  JDATA  = 0^  DTAU 
is  automatically  calculated  within  the  program. 

Computation  to  proceed  in  time  up  to  this  value. 

Nominal  time  level  at  start  of  program.  Output  is  indexed 
with  respect  to  this^  so  as  to  have  output  at  convenient 
time  values. 

relaxation  parameter  for  iteration  of  stream  function 
field . 

If  = 0^  calculation  starts  from  quiescent  initial 
conditions'  if  it  starts  from  input  data  fields. 

maximum  number  of  stream  function  iteration  cycles 
permitted  at  each  time  step.  Each  cycle  consists  of 
two  iterative  sweeps  over  the  field;  one  vertically 
and  one  radially. 

heat  transfer  rates  calculated  and  printed  at  times  TAU  - 
TAUNOM  + n * FKTR^  where  n =1  0^  1^  2^  • * • . 

fields  of  temperature^  velocity^  vorticity  and  stream 
function  printed  at  times  TAU  = TAUNOM  + n • FPRINT_, 
v;here  n = 0^  1^  2^  * * • . 

fields  put  on  magnetic  tape  at  times  TAU  = TAUNOM  + 
n • FTAPE_,  where  n = 0^  1^  2^  • * * . 

fields  punched  onto  cards  at  times  TAU  = TAUNOM  + n * FPUNCH^ 
where  n = 0^  1^  2^  • • • . 
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= not  used  in  present  program. 

= logical  I/O  unit  on  which  input  tape  is  mounted  (if  used) . 

= logical  I/O  unit  on  which  output  tape  is  mounted  (if  used). 
= data  block  to  be  read  from  input  tape. 

= data  block  in  which  writing  starts  on  output  tape. 


After  reading  the  input  parameters  the  initial  data  fields  are  either 
generated  within  the  program  or  read  in^  according  to  the  value  of  JDATA. 
The  program  then  enters  the  computation  cycle  as  shown  in  figure  2. 


2.2  Subroutines 


SUBROUTINE  VSL  (p  8 of  figure  3):  This  subroutine  is  called  from  the 

main  program  and  computes  the  velocities  U and  V at  node  points  with  central 
difference  approximations.  At  the  solid  boundaries^  U = V = 0.  The  output 
fields,  however,  have  nonzero  values  for  U along  the  vertical  wall  and  for  V 
along  the  floor  and  ceiling.  These  values  pertain  to  a location  AR/2  or 
Ax/2  away  from  the  boundary  and  are  computed  with  one  sided  differences. 

SUBROUTINE  MEAN  (p  9 of  figure  3):  This  subroutine  is  called  from 

the  main  program  and  computes  the  following  "mean”  velocities: 

UM(I,J)  = (U(I,J)  + U(I+1,  J))/(2AX) 

VM(I,J)  = (V(I,J)  +V(I,  J+1))/(2AR). 

The  grid  spacings  AX  and  Ar  appear  in  the  denominators  and  eliminate 
several  multiplications  in  the  main  program. 

SUBROUTINE  RESETS  (p  10  of  figure  3):  This  subroutine  is  called  from 

the  main  program  and  controls  printing  of  the  fields  of  temperature,  velocity, 
vorticity  and  stream  function. 

SUBROUTINE  BLOCK  (p  11  of  figure  3):  This  is  called  from  subroutine 

RESETS  to  break  up  and  print  the  output  fields  in  blocks  ten  or  eleven 
columns  wide.  The  top  and  bottom  of  each  block  corresponds  to  the  ceiling 
and  floor,  respectively.  The  printed  output  is  illustrated  in  the  example, 
section  3. 

SUBROiri'INE  PUN  (p  12  of  figure  3):  This  subroutine  is  called  from  the 

main  program  and  punches  the  fields  of  temperature,  velocity,  etc.,  onto 
cards.  The  first  card  serves  to  identify  the  data  field.  When  used  as  input 
for  the  main  program,  this  first  card  replaces  the  first  of  the  three  cards 
(shown  on  p 16  of  figure  3)  which  contain  the  input  parameter  list.  This 
first  card  is  then  followed  by  two  cards  containing  suitable  values  for  the 
remaining  input  parameters,  and  the  data  field  cards. 
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SUBROUTINE  TAP  (p  13  of  figure  3):  When  called  from  the  main  program_, 

this  subroutine  puts  the  fields  of  temperature,  velocity,  etc.,  onto 
magnetic  tape,  followed  by  an  end  of  file  mark.  (An  end  of  file  mark  is 
always  placed  at  the  beginning  of  a new  tape).  The  first  record  identifies 
the  data  block  and  is  of  the  same  format  as  the  first  input  data  card  of 
the  main  program.  This  record  is  used  to  locate  data  blocks  on  tape.  The 
information  in  this  record,  in  the  same  format,  is  also  punched  onto  two 
cards . 

SUBROUTINE  CHECK  (p  14  of  figure  3):  VThen  called  from  the  main  program, 

this  subroutine  scans  the  first  record  after  every  end  of  file  mark  on  the 
input  tape,  until  a data  block  corresponding  to  the  first  input  data  card 
of  the  main  program  is  located.  Four  input  data  cards  are  required  for  the 
main  program,  the  first  and  fourth  cards  being  identical. 

SUBROUTINE  NUSSEL  (pp  15-16  of  figure  3):  When  called  from  the  main 

program,  this  subroutine  calculates  two  rates  of  heat  transfer:  (a)  from 

heat  source  to  fluid  in  enclosure;  and  (b)  from  fluid  in  enclosure  to 

cold  walls.  Results  are  printed  as  a dimensionless  ratio,  /XaAT, 

where  \ is  the  thermal  conductivity  of  the  fluid  in  the  enclosure. °^^uss el t 

numbers  are  also  computed  by  dividing  the  dimensionless  rates  of  heat 

transfer  by  the  dimensionless  heating  area  (A,  /a*')  or  cooling  area  (A  /a  ). 

in  out 

3 . EXAMPLE 


Results  obtained  from  the  execution  of  the  program  sho\>m  in  figure  3 
are  presented  in  figure  4.  Page  1 of  the  figure  is  the  heading  for  the 
output  results,  which  lists  the  pertinent  input  parameters  for  the  problem. 
These  are  some  of  the  parameters  read  in  from  the  data  cards  on  p 16  of 
figure  3,  Following  this  hdading,  heat  transfer  results  and  calculated 
fields  are  periodically  printed  out.  To  save  space,  only  the  printed  out- 
put for  TAU  = 0.25  is  shown  (pp  2-6  of  figure  4).  The  computed  fields  are 
listed  in  the  following  order:  temperature,  vertical  velocity  U,  radial 
velocity  V,  vorticity,  and  stream  function.  A 21x21  grid  V7as  employed, 
and  for  each  field,  441  numbers  are  printed.  For  the  first  block  of 
numbers,  21x10,  the  left-most  column  pertains  to  the  centerline,  and 
columns  to  the  right  correspond  to  increasing  values  of  radius,  up  to 
R = 0,45.  The  top  and  bottom  numbers  in  each  column  correspond  respectively 
to  the  ceiling  and  the  floor.  The  second  block  of  numbers,  21x11,  covers 
the  radius  range  from  R = 0,5  (left  column)  to  R = 1.0  (right  column). 

Again,  ceiling  and  floor  correspond  to. the  top  and  bottom  of  each  column. 

As  noted  in  the  discussion  of  SUBROUTINE  VEL,  the  right-most  column 
of  vertical  velocities  U is  nonzero.  These  values  pertain  to  R = 0.975 
(not  R = 1,0).  Similarly,  the  radial  velocities  V along  the  floor  and 
ceiling  are  nonzero.  These  values  respectively  correspond  to  X = 0.025 
and  X = 0.975. 
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Following  the  stream  function  field_,  some  additional  output  is  shown 
(p  6 of  figure  4).  The  first  entry  is  the  number  of  cycles  required  for 
convergence  of  the  stream  function  iteration  (in  this  case_,  one).  The 
next  three  entries  constitute  printed  output  from  the  subroutines  PUN,  TAP, 
and  NUSSEL.  Similar  output  results  results  whenever  these  subroutines  are 
called.  The  words  CARDS  PUNCHED  indicate  that  the  temperature,  velocity, 
etc.,  fields  have  been  punched  onto  cards.  The  next  line  indicates  that 
these  fields  have  been  written  on  tape  in  BLOCK  9.  The  first  identifying 
record  of  the  data  block  on  tape  is  printed  between  the  asterisks.  The 
remaining  two  lines  on  the  page  are  printed  by  subroutine  NUSSEL,  and 
provide  the  time  TAU,  the  next  time  step  DTAU,  and  the  heat  transfer 
rates . 

4.  APPLICATION  EXPERIENCE 


Most  of  the  calculations  to  date  with  EXCONl  have  employed  fixed 
values  of  the  following  parameters:  RB  = 1,  RS  = 0.1,  and  PR  = 0.7.  This 

includes  reference  2,  which  presents  information  on  ^runcation  ^rrors  and 
the  effects  of  grid  size.  A range  of  GR  from  1 x 10^  to  4 x 10^  is  ^overed. 
Referejge  3 presents  a study  of  the  effect  of  varying  GR  from  4 x 10  to 
4 X 10  with  M = 50  and  N = 20.  The  computed  flows  are  all  laminar.  The 
validity  of  the  calculations  for  GR  < 1 x 10  is  confirm.ed  by  comparison 
with  the  experiments  of  reference  4.  At  higher  values  of  GR  the  experiments 
indicate  turbulence,  although  the  calculated  flows  are  laminar.  This 
difference  is  due  to  the  assumed  axisymmetry  and  to  numerical  truncation 
errors.  The  latter  increase  with  both  GR  and  coarser  grids.  From  figure 
5,  it  is  apparent  that  the  M x N = 50  x 20  grid  employed  in  reference  3 
required  computing  times  of  several  hours  at  high  values  of  GR. 

g 

Some  calculations  have  also  been  attempted  at  GR  = 4 x 10  with  an 
aspect  ratio  of  RB  = 10.  Preliminary  runs  with  10  x 100  and  20  x 200  grids 
indicate  that  the  required  computing  times  to  reach  steady  state  are  81 
minutes  and  17  x 81  minutes,  respectively.  Both  grids  are  too  coarse  to 
adequately  resolve  details  of  the  flow.  Furthermore,  the  finer  grid 
increases  the  amount  of  computer  time  required  roughly  as  the  ratio  of  the 
squares  of  the  number  of  grid  points. 

From  experience,  the  following  conclusion  can  be  noted  about  the  use 
of  EXCONl:  Meaningful  results  can  be  obtained  only  when  the  grid  is  fine 

enough  to  resolve  details  of  the  flows.  If  the  grid  is  too  coarse,  details 
are  not  resolved  and  truncation  errors  become  important.  The  use  of  fine 
grids  to  improve  the  calculation,  however,  can  quickly  lead  to  prohibitively 
large  demands  on  com.puter  time  and  storage.  The  application  of  EXCONl  to 
very  high  GR  natural  convection  flows  does  not  yet  seem  feasible,  because 
of  the  foregoing  requirements.  Nevertheless,  the  program  can  be  effectively 
and  successfully  applied  to  calculate  flows  when  the  mesh  size  is  chosen 
judiciously  and  does  not  exceed  computer  storage  limitations. 
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APPENDIX  A.  APPLICATION  TO  THE  FACTORY  ^!UTUAL 
IBM  360  COMPUTER  IN  NORWOOD,  MASSACHUSETTS 

The  program  listed  in  figure  3 must  be  modified  for  use  on  the  FM 
IBM  360.  The  IBM  compiler  uses  a ’’basic  Fortran”  language  and  does  not 
recognize  logical  IF  statements^  which  are  used  extensively  throughout 
EXCONl.  These  statements  would  all  have  to  be  reprogrammed.  An  improved 
IBM  compiler  is  being  developed^  which^  when  available^  will  recognize 
the  logical  IF.  Only  minor  modifications  would  then  be  required. 

Any  decks  of  cards  punched  for  or  on  the  Univac  1108  have  to  be 
converted  for  use  on  the  IBM  360.  A fev7  symbols  are  punched  differently 
for  the  two  computers^  such  as  parentheses^  plus  sign^  etc.  Dr.  John  A. 
Rockett  has  developed  a program  to  effect  such  a conversion. 
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FIGURE  1.  Cylindrical  room  and  grid  system. 


FIGURE  2.  Flowchart  of  natural  convection  calculation.  The  corresponding 
computer  program  is  shown  in  figure  3. 


FIGURE  3 (16  pages).  Listing  of  EXCONl 
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Fortran  language  program. 
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yiOVRE  -t  (a  pages),  Reiu-csoi'tative  printed 


outpvit  Cor  example  case 


\iT0i?4L  CONVECTION  IN  A ROOM  ARISING  FROM  A LOCALIZED  KEAT  SOURCE  AT  FLOOR  LEvEL- 
HEAT  SOURCE  AT  TEMPERATURE  T . 1, 

CEILING,  A alls,  ANO  remainder  OF  FLOOR  AT  T = 0, 

INPUT  parameters 

RADIUS/HEIGhT  3 B/A  = RB  = I.OqO  HEAT  SOURCE  R A 0 1 US / R A D I US  = C/B  i RS  = ,iQO 

GRASHOF  BASED  ON  HEIGHT  = gR  = 1.000+05  PRANDTL  = PR  = .700 

M = 20  N = 20  0*  = .050000  DR  * .050000 

TAUZER  = 0.0000000  DTaUZE>^0  = 3.9583333-OR  TAUMAX  = .250000QO 

RELAX  = 1.75  ITMAX  = 15  kTr  • .0010000  PRINT  = .0100000  TAPE  = .0100000  PUNCH  = .2500 


1 nt 


I'lGURR  4 (concM) 


TAU  ■ 2.506*1057-0  1 


THE 

000000 

TEMPERATURE 

.000000 

field  Is 
. DOOOOO 

G I VEN  BY 
.000000 

.000000 

.000000 

,000000 

.000000 

. qoqooo 

.000000 

0 1 8008 

. 0 1 7392 

.016278 

•013639 

.012827 

.01 1022 

.009331 

.007809 

.006377 

.005338 

031363 

.030607 

.028581 

•025868 

.022882 

.019901 

.017080 

.013396 

.012177 

• 0 1 0 1 2 3 

030075 

.038952 

.036238 

•032699 

. 028906 

.025229 

.021828 

.018738 

.015976 

• 0 1 3 3 8 3 

035529 

.033157 

.030686 

•036233 

.031657 

.027325 

.023693 

.020339 

.017603 

• 0 1 5067 

039712 

.037919 

•033513 

•037979 

.032396 

.027667 

.023633 

•020336 

.017573 

• 0 1 5 1 85 

053551 

.051306 

•035786 

•038983 

.032363 

.026962 

.022616 

.019236 

•016588 

• 0 1 3 39  7 

057682 

.053890 

.038058 

• 0398  1 1 

.032137 

.025926 

.021232 

. 13  1 7 7 7 8 

.015200 

• 0 1 3 1 83 

062360 

.058998 

.050610 

•030723 

.031833 

.023867 

.019825 

.016275 

.013736 

• 0 1 1 86 1 

068183 

.063880 

.053603 

• 03 1 83  1 

.031633 

.023953 

.018572 

.013905 

•012388 

•010590 

0751 99 

.069797 

.057210 

• 033307 

.031691 

.023239 

.017398 

.013699 

.011172 

• 009328 

083998 

.077163 

• 06  1 682 

•035220 

.03i99i 

.022700 

.016580 

.012633 

.010093 

.008387 

095277 

.0865  1 1 

.067253 

• 03762  1 

.032510 

.022293 

.015779 

.01 1702 

.009  1 29 

•007352 

110060 

.098606 

.073259 

• 050580 

.033220 

.021973 

.015037 

.010837 

.008237 

•006600 

129932 

.113676 

.083230 

•053206 

.033092 

.021679 

.013331 

.010003 
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FIELD  OF  VELOCITY  COMPONENT  U 

IS  GIVEN  BY 

0.000 
5.112+00 
I . M54  + 0 1 
2.335+01 
3.  139  + 01 


0.000 
3.733+00 
1 .363*01 
2.235*01 
2.963*0 1 


0.000 
3.086*00 
1 . 202*01 
2.005*01 
2 .659*01 


0.000 
3.220*00 
9.832*00 
1 .673*01 
2.232*01 


0.000 
2 . 3 I 9*00 
7 .383*00 
1 . 3 1 0*0 1 
1 .783*01 


0 .000 
1.508*00 
5.219*00 
9.675*00 
1 .330*01 


0.000 
8.505-01 
3 .276600 
6.327+00 
9 .393+00 


0.000 
3.607-01 
1 .710*00 

3.773*00 

5.936*00 


0.000 

2.333-02 

5.298-01 
1 .639*00 
3.065*00 


0.000 
-1  .866-0  1 
-2.953-01 
3.915-02 
7.683-01 


|G£! 


3,494.:1 
M , 092 ♦!  1 
M.332-.C1 
9 . M 7 7 • 1 

m.5‘U»31 
M . S9  3 *0  1 
>1 . 4 96  *0  1 
4 . 404 ♦ 0 1 
4.269*01 
4, 083*0  1 
3.829*01 
3.471*01 
2.955*01 
2.199*01 
1 .096*01 
0 .000 

0.000 
-3.096-01 
-8.292-01 
-1.1 00*00 
-9.752-01 
-5.716-01 
-9 . 024-02 
3. 138-0 1 
5.634-01 
6.437-01 
5.788-0 1 
4. 104-01 
1.851-01 
-5.342-02 
-2.659-01 
-4  . 1 49-0 1 
-4.780-01 
-4.323-01 
-2.884-01 
-1 .021-01 
0.000 


THE 


0.000 
0.000 
0 .000 
0.000 
0.000 
0.000 
0.000 
0 .000 
0.000 
0.000 
0.000 
0.000 
O.QOO 
O.QOO 
0.000 


3.490*01 
3.848*01 
4.070*01 
4 . 188*01 
4.223*01 
4 . 1 95*0 1 
4.114*01 
3.991*01 
3.820*01 
3.597*01 
3 . 309*01 
2.929*01 
2 .420*0 1 
1.731*01 
8 .340*00 
0.000 


3.134*01 
3.446*01 
3.426*01 
3.702*01 
3.698*01 
3.630*01 
3.51 1*01 
3.346*01 
3. 135*01 
2.875*01 
2.556*01 
2.164*01 
1 .684*01 
1 . 1 05*01 
4 .622*00 
0.000 


2.449*01 
2.904*01 
3.036*01 
3 .069*01 
3.024*01 
2.924*01 
2.773*01 
2.580*01 
2 .344*01 
2.071*01 
1.751*01 
1 .383*01 
9 .498*00 
5 . 3 16*00 
1 .438*00 
0.000 


2. 123*01 
2.328*01 
2.420*01 
2.424*01 
2 . 359*0 1 
2.242*01 
2.082*01 
1 .987*01 
1 .659*01 
1 .402*01 
1 . 1 1 8*0 1 
8. l05*Oo 
4 .945*00 
2.028*00 
4.454-02 
0 .000 


1 . 6 18*0 1 

1 .784*01 

1 .853*01 

1 .844*01 

1.774*01 
1.642*01 
1 .513*01 
1 .334*0 1 
1 . 136*0 1 
9. 187*00 
6 .886*00 
4.556*00 
2.344*00 
6.312-01 
-1 .989-01 
0.000 


1.144*01 
1 . 302.0  1 
1 . 360*01 
1 .352*01 
1 .293*01 
1 . 194*0 1 
1 .070*01 
9.223*00 
7 .599*00 
5.883*00 
4. 147*00 
2.491 *00 
1.072*00 
1 .087-01 
-1 .882-01 
0.000 


7 . 708*00 
8.872*00 
9.415*00 
9 . 420*00 
8.996*00 
8.248*00 
7 .247*00 
6.126*00 
4.891 *00 
3.622*00 
2.390*00 
1 .283*00 
4 . 142-01 
.7.943-02 
-1 .427-01 
0.000 


4.379*00 

5 .342*00 
5.849*00 
5.979*00 
5.743*00 
5.240*00 
4.549*00 
3.739*00 
2.871*00 
2.003*00 
1 . 197*00 
5.196-01 
4. 408-02 
-1 .404-01 
-1.113-01 
0.000 


1.  627*00  FIOVRli 
2 . 368*00 
2.857*00 
3.058*00 
2.997*00 
2 .728*00 
2.313*00 
1 .81 1*00 
I .276*00 
7.598-01 
3.093-01 
-2.821-02 
-2. 1 10-01 
-2.203-01 
-9.879-02 
0.000 
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-4.032-0 1 
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-2.2 15*00 
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-2.240*00 

-3.428*00 

-2.050*00 

-3.531*00 

- 1 . 785*00 

-3.446*00 

-1 .545*00 

-3.271*00 

- 1 , 374*00 

-3.045*00 

- 1 . 277*00 

-2.858*00 

-1 .237*00 

-2.653*00 

- 1 . 222*00 

-2.440*00 

- 1 . 204*00 

-2.204*00 

-1.1 54*00 

-1.931 *00 

-1.05 1*00 

-1.41 3*00 

-8.810-01 

- 1 , 248*00 

-6.456-01 

-8.5 16-0 1 

-3.497-01 

-4 . 585-0 1 

- 1 . 168-0  1 

-1,379-01 

0.000 

0.000 

0.000 

0.000 

-4.077-01 

-3.992-01 

- 1 . 380*00 

- 1 . 380*00 

-2.522*00 

.2.597*00 

-3.516*00 

-3.749*00 

-4 . 204*00 

-4.498*00 

-4.575*00 

-5.231*00 

-4.663*00 

-5,514*00 

-4.609*00 

-5.554*00 

-4.416*00 

-5 . 4O 1 *00 

-4 . 1 44*00 

-5. 107*00 

-3.822*00 

-4 , 706*0q 

-3. 454*00 

-4,224*00 

-3.041*00 

-3.677*00 

-2,585*00 

-3.082*00 

-2.090*00 

-2.452*00 

-1 .566*00 

-1 ,807*00 

-1 .036*00 

- 1 . 177*00 

.5.433-01 

-6.099-01 

- 1 . 405-0 1 

-1 .793-01 

0.000 

0.000 

0.000  0.000 
-3.833-01  -3.592-01 
-1.338*00  -1.255*00 
-2.543*00  -2.430*00 
-3.784*00  -3.642*00 
-4.808*00  -4.704*00 
-5,539*00  -5.503*00 
-5.958*00  -5,997*00 
-4.093*00  -6, 192*00 
-5.988*00  -6.122*00 
-5.691*00  -5.834*00 
-5,244*00  -5.375*00 
-4,692*00  -4,791*00 
-4.059*00  -*♦ . 122*00 
-3.372*00  -3.399*00 
-2.655*00  -2.655*00 
-1.9l7*00  -1.919*00 
-1,250*00  -1,227+00 
-6.425-01  -4.240-01 
-1 .887-01  -1 .834-01 
0.000  0.000 


O.QOO  0.000 

-3,174-01  -2.405-01 
-1.111 *00  -8.635-01 
-2,170*00  -1.71 9*00 
-3.290*00  -2.642*00 
-4.297*00  -3.482*00 
-5.075*00  -4,134*00 
-5,570*00  -4,547*00 
-5.777*00  -4.713*00 
-5,721*00  -4,656*00 
-5,448*00  -4.415*00 
-5,005*00  -4,034*00 
-4,442*00  -3.559*00 
-3,798*00  -3,022*00 
-3, 1 10*00  -2.454*00 
-2.408*00  -1,880*00 
-1.722*00  -1.324*00 
-1.088*00  -8,212-01 
-5,473-01  -4.010-01 
-1.583-01  -1.105-01 
0,000  0.000 


0.000  0.000 
-1.159-01  -4.293-02 
-4.741-01  -2.443-01 
-9.917-01  -5.540-01 
-1 .558*00  -8.958-01 
-2.073*00  -1.204*00 
-2.467*00  -1.434*00 
-2.708*00  -1.569*00 
-2.795*00  -1,409*00 
-2.745*00  -1.569*00 
-2.584*00  -1.468*00 
-2.347*00  -1.323*00 
-2.055*00  -1.150*00 
-1 .732*00  -9.619-01 
-1.394*00  -7,479-01 
-1.054*00  -5.750-01 
-7.315-01  -3.904-01 
-4.384-01  -2.242-01 
-2.002-01  -9,122-02 
-4,622-02  -1 .231-02 
0,000  0.000 


FIELD  OF  velocity  COMPONENT  V IS  GIVEN  BY 


2 . 556*00 
3 .440*00 
4 , 485*00 
4 .207*00 
3 .277*00 
2.358*00 
1 . 590*00 
9.863-01 
5 ,220-0  I 
1 . 658-0 1 
- 1 . 140-0 1 
-3.475-01 
-5,451-01 
-8.020-01 
-1 . 102*00 
-1 .530*00 


4 ,733*00 
4 . 8 13*00 
8 . 860*00 
8 .002*00 
6 .221*00 
4 ,424*00 
2 .904*00 
1 . 699*00 
7,444-01 
3.935-02 

-5,374-01 
-1 ,024*00 
-1.481 *00 
- 1 . 969*00 
-2.540*00 
-3. 344*00 


6 . 300*00 
9.228*00 
1 .220*01 
1 . 1 12*01 
8.621*00 
6.032*00 
3.81 1*00 
2 .036*00 
6.536-0 1 
-4.234-01 
- 1 . 283*00 
-2.01 3*00 
-2.695*00 
-3 . 408*00 
-4.230*00 
-5.246*00 


7 , 197*00 
1 .078*01 
1 , 456*0 1 
1 .384*01 
1.043*01 

7 . 1 84*00 
4 . 355*00 
2.081*00 
3.085-01 

-1,068*00 

-2 . 1 65*00 
-3 . Q90*00 
-3.941*00 
-4.801 *00 
-5.745*00 

-6,830*00 


7,490*00 

1.151*01 

1 .594*01 

1.497*01 
1 . 148*01 

7 , 970*00 
4 , 645*00 
1 .987*00 
-9 , 880-02 
- 1 .7  12*00 
-2.987*00 
-M, 048*00 
-«. 998*00 
-5.919*00 
-6.871 *00 
-7,881 *00 


7.311*00 
1 . 154+01 
1 .643*01 
1 ,578*01 
1 ,244*01 
8,489*00 
4 .860*00 
1 ,885*00 
-4.387-01 
-2.230*00 
-3.632*00 
, 776*00 
-5 ,746*00 
-6.679*00 

.7,542*00 

-8,414*00 


4 ,808*00 
1.103*01 
1 .617*01 

1 ,594*01 
1 ,281*01 
8.804*00 
5,012*00 
1 .847*00 
.6,436-01 
-2,562*00 

-4.050*00 

-5,237*00 

-4,224+00 

-7.091 *00 
.7,867*00 
-8.534*00 


4,114*00 
1 .015*01 
1 .531*01 
1 .553*01 
1.279*01 
8,936*00 
5. 139*00 
1 ,892*00 
-6.961-01 

-2.497*00 

-4.237*00 
-5 . 439*00 
-4.404*00 
-7.201*00 
-7 ,859*00 
-8 . 357*00 


5.339*00 
9.048*00 
1 .402*01 
1 .445*01 
1 .239*01 
8.872*00 
5.222*00 
2 .003*00 
.4. 125-01 
-2.450*00 
-4.21 1*00 
-5,408*00 
-6.334*00 
-7.058*00 
-7 .406*00 
-7 ,958*00 


o.oijo 

0 t 000 
0.000 
0.000 
0.000 

4 . 554  »0Q 
7 .853*00 
1 .245*01 

1 . 33a*0 1 
1 . I 64*01 
8.589*00 
5.218*00 
2.  1 35*00 

-4.305-01 
-2.455*00 
-4.005*00 
-5. 177*00 
-6,056*00 
-4.706*00 
-7 , 1 58+00 
-7.395*00 
-7.341*00 
-6.840*00 
-5.436+00 
-3,387*00 
-1 .954*00 


THE 


-2 . 185*00 
-3.204*00 
-4,446*00 
-5.473*00 
-5.431*00 

3 .805*00 
4 .449*00 
I .075*01 
1.185*01 
1 .042*01 
8 .063*00 
5.077*00 
2.233*00 
-2 . 0 12-0  I 
-2. 154*00 
-3.657*00 
-4,782*00 
-5.605*00 
-4.1 84*00 
-6.553*00 
-4.704*00 
-4.578*00 
-4.049*00 
-4.920*00 
-2 .925*00 
- 1 . 475*00 


F I ELO  OF 


-4.440*00 
-5.985*00 
-7.928*00 
-3.455*00 
-8 . 340*00 

3 . 103*00 
5 .492*00 
9.024*00 
1.018*01 
9 .355*00 
7 .308*00 
4 .765*00 
2 .242*00 
1 .850-02 
- 1 . 800*00 
-3 . 2 12*00 
-4.264*00 
-5.017*00 
-5,526*00 
-5.824*00 
-5 . 9 12*00 
-5.747*00 
-5.231*00 
-4 , 2 12*00 
-2.484*00 
-1 .414*00 


-4.538*00 
-8 . 130*00 
-9.496*00 
-9 . 192*00 
.7 . 990*00 

2.475*00 
4.412*00 
7 .352*00 
8.445*00 
7.940*00 
4.342*00 
4 .279*00 
2.129*00 
1 . 8 19-01 

- 1 . 44  2*00 
-2.714*00 
.3.441 *00 
-4 . 328*00 
-4 . 743*00 
-5.000*00 
-5.043*00 
.4.843*00 
-4.389*00 
-3 . 504*00 
-2.052*00 

- 1 . 162*00 


S I VEN 


-8 .079*00 
-9.379*00 
-1 .014*01 
-8 . 3 15*00 
-4.327*00 

1 .907*00 
3 . 424*00 
5.774*00 
4.738*00 
4.455*00 
5 .282*00 
3 .641*00 
1 .890*00 
2 .642-01 
-1.1 13*00 
-2.202*00 
-3.01 1*00 
-3,574*00 
-3.930*00 
-4 . 109*00 
-4.121 *00 

-3.947*00 

-3.534*00 

-2.797*00 

- 1 .425*00 

-9 . 142-0 1 


-8 .908*00 
-9.740*00 
-9. 720*00 

-7 . 138*00 

-4 ,9oO*00 

1 . 400*00 
2.536*00 
4 ,320*00 
5. 109*00 
4 . 948*00 
4 . 1 30*00 
2.898*00 
1.545*00 
2.644-01 
-8.323-01 

- 1 . 703*00 
-2.348*00 
-2.791*00 
-3.043*00 
-3.189*00 
-3.1 80*00 
-3.025*00 
-2.687*00 
-2. 105*00 

- 1 . 209*00 
-6.727-01 


0,000 

-1.3 17*02 

-2 

.424*02 

.3 

. 195*02 

-3,602*02 

-3.490*02 

0.000 

-3.328*0 1 

-6 

,450*01 

-9 

. 8 12*01 

- 1 , 253*02 

-1  .454*02 

0.000 

2.455*01 

.260*01 

4 

.449*01 

4,030*01 

2.454*01 

0.000 

5 .488*01 

9 

.521*01 

1 

. 1 90*02 

1 .275*02 

1 .241*02 

0.000 

4.442*01 

1 

. 1 66*02 

1 

.473*02 

1 .414*02 

1 .430*02 

0.000 

7.218*01 

1 

.254*02 

1 

.544*02 

1 , 494*02 

1 .705*02 

O.QOO 

7.421*01 

1 

.306*02 

1 

.590*02 

1 .674*02 

1 .449*02 

0.000 

8.037*0  1 

1 

. 354*02 

1 

. 40 1 *02 

1 .632*02 

1 .555*02 

0.000 

8,521*01 

1 

, 407*02 

1 

. 4 1 4*02 

1 .586*02 

1 .454*02 

O.QOO 

9.095*01 

1 

.469*02 

1 

.631*02 

1 .544*02 

1 .342*02 

0.000 

9 .775*01 

1 

.540*02 

1 

.654*02 

1 .51 1*02 

1 .290*02 

0.000 

1 .054*02 

1 

.621*02 

1 

.684*02 

1 . 484*02 

1 .222*02 

0 .000 

1 . 147*02 

1 

. 7 13*02 

1 

. 7 19*02 

1 .440*02 

1 . 157*02 

0,000 

1 ,251*02 

1 

. 8 15*02 

1 

.755*02 

1.43 1*02 

1 .089*02 

0.000 

1 .367*02 

1 

.923*02 

1 

.784*02 

1 .390*02 

1.008*02 
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